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Abstract: We perform the Hamiltonian analysis of non-relativistic covariant Hofava- 
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relativistic covariant Hofava-Lifshitz gravity presented in arXiv:1007.2410. Then we extend 
this construction to the case of RFDiff invariant Hofava-Lifshitz theory. We find well 
behaved Hamiltonian system with the number of the first and the second class constraints 
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1. Introduction and Summary 

In 2009 Petr Hofava formulated new proposal of quantum theory of gravity that is power 
counting renormalizable |3|. This theory is now known as Hofava-Lifshitz gravity 

(HL gravity). It was also expected that this theory reduces do General Relativity in the 
infrared (IR) limit. HL theory was studied from different point of view due to the fact 
that this is a new and intriguing formulation of gravity as a theory with reduced amount 
of symmetries that leads to remarkable new phenomena . 

The HL gravity is based on an idea that the Lorentz symmetry is restored in IR limit 
of given theory and can be absent at high energy regime of given theory. Explicitly, Hofava 
considered systems whose scaling at short distances exhibits a strong anisotropy between 
space and time, 

x' = be , t' = Ft . (1.1) 

In (D + 1) dimensional space-time in order to have power counting renormalizable theory 
requires that z > D. It turns out however that the symmetry group of given theory 
is reduced from the full diffeomorphism invariance of General Relativity to the foliation 
preserving diffeomorphism 

x" = x i + ( i {t,x.) , t' = t + f(t). (1.2) 

Due to the fact that the diffeomorphism is restricted ( |1.2|) one more degree of freedom 
appears that is a spin— graviton. It turns out that the existence of this mode could 
be dangerous since it has to decouple in the IR regime, in order to be consistent with 
observations. Unfortunately, it seems that this might not be the case. It was shown that 
the spin-0 mode is not stable in the original version of the HL theory [1] as well as in the 
Sotiriou, Visser and Weinfurtner (SVW) generalization ||. Note that in both of these two 
versions, it was all assumed the projectability condition that means that the lapse function 
N depends on t only. This presumption has a fundamental consequence for the formulation 
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of the theory since there is no local form of the Hamiltonian constraint but the only global 
one. However we would like to stress that these instabilities are all found in the Minkowski 
background. Recently, it was found that the de Sitter spacetime is stable in the SVW setup 
H [l(]]. Then we can presume that this background is legitimate background. 

On the other hand there is the second version of HL gravity where the projectability 
condition is not imposed so that N = N(x,t). Properties of given theory were extensively 
studied in 0, |l2|, |l|, [U|, [H], |l|, ^ |(| |l], . It was shown recently in @ that so 

called healthy extended version of given theory could really be an interesting candidate for 
the quantum theory of reality without ghosts and without strong coupling problem despite 



its unusual Hamiltonian structure [17, 18] 



Recently Hofava and Malby-Thompson in [25| proposed very interesting way how to 
eliminate the spin-0 graviton. They considered the projectable version of HL gravity to- 
gether with extension of the foliation preserving diffeomorphism to include a local U(l) 
symmetry. The resulting theory is then called as non-relativistic covariant theory of grav- 
ity 2 . It was argued there |H|] that the presence of this new symmetry forces the coupling 
constant A to be equal to one, however this result was questioned in ]^6| (see also fl3C| l) 
where an alternative formulation of non-relativistic general covariant theory of gravity was 
presented. Further, it was shown in [25^, 26] that the presence of this new symmetry implies 
that the spin-0 graviton becomes non-propagating and the spectrum of the linear fluctua- 
tions around the background solution coincides with the fluctuation spectrum of General 
Relativity. 

This new proposal of non-relativistic general covariant HL gravity is very interesting 
and it certainly deserves further study. In this paper we present the Hamiltonian analysis 



of the formulation of non-relativistic covariant HL gravity given in [26|. We argue that re- 



sulting Hamiltonian and constraint structure has the same form as in [25] even if they differ 
in explicit form since they are derived from different Lagrangians. This fact shows that 
these two Lagrangian formulations of non-relativistic covariant HL gravities are equivalent 
on the level of the Hamiltonian formalism as well. 

Despite the fact that non-relativistic covariant HL gravity seems to solve the problem 
of the scalar graviton and the content of the physical degrees of freedom is the same as 
in General Relativity there is still one additional first class constraint which is the global 
Hamiltonian constraint. The meaning of this constraint should be investigated further as 
was nicely discussed on page 30 in 2E]. In order to find version of non-relativistic covariant 
HL gravity without global Hamiltonian constraint we recall that there exists formulation 
of the HL gravity with reduced symmetry group known as restricted-foliation-preserving 
Diff (RFDiff) HL gravity 14, 23]. This is the theory that is invariant under following 
symmetries 



t' = t + St ,5t = const 



x* + C < (x ) t) • 



(1-3) 



The characteristic property of given theory is that in its simplest version [^] based on the 
detailed balance construction M, 0, [| there is no reason to introduce the lapse function 



2 This theory was also studied in ||| fj|, ^ Q. 
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N 3 . Then we introduce U(l) symmetry as in |2f| or its alternative version given in |p6| . 
Finally we proceed to the Hamiltonian formulation of given theory and we find there is no 
global Hamiltonian constraint due to the absence of the lapse function N in the action. We 
further determine all constraints in given theory and we show that the number of the first 
class and the second class constraints implies that the physical phase space has the same 
dimensions as in case of General Relativity. On the other hand we show that the presence 
of the second class constraints implies that the symplectic structure of given theory that is 
determined by corresponding Dirac brackets between physical degrees of freedom is rather 
complicated due to the fact that Dirac brackets generally depend on phase space variables. 

Let us outline our results and suggest possible extension of this work. We perform 
the Hamiltonian analysis of the theory suggested in [26| and we show that its Hamiltonian 
structure is equivalent to the Hamiltonian structure found in paper p5[|. We also suggest 
an alternative formulation of non-relativistic covariant HL gravity that is based RFDiff HL 
gravity. We show that the resulting theory has consistent Hamiltonian formulation with 
the same content of the local constraints as in case of non-relativistic covariant HL gravity 
but without global Hamiltonian constraint. On the other hand we should stress that we 
are not able to solve explicitly the second class constraints with respect to physical degrees 
of freedom in the full generality. We are also not able to determine corresponding Dirac 
brackets. Then it would be clearly desirable to find exact results at least for some special 
situations. It would be also interesting to find exact solutions of the equations of motion of 
non-relativistic general covariant RFDiff HL gravity. We hope to return to these problems 
in near future. 

The organization of this paper is as follows. In the next section (||) we introduce the 
non-relativistic general covariant HL gravity in the formulation firstly presented in [26]. 
Then in section @ we perform its Hamiltonian analysis. In section (||) we introduce the 
non-relativistic general covariant RFDiff-invariant HL gravity. Then we perform its Hamil- 
tonian analysis and shows that the resulting theory correctly describes physical degrees of 
freedom of D + 1 dimensional gravity. 



2. Non-Relativistic Covariant HL Gravity 

We begin this section with the introduction of basic notation, for detailed treatment of 
D + 1 formalism, see [24]] . 

Let us consider D + 1 dimensional manifold M with the coordinates , fj, = 0, . . . , D 
and where x^ = (t,x) , x = (x 1 , . . . , x D ). We presume that this space-time is endowed 
with the metric g^ u (x p ) with signature (—,+,...,+). Suppose that M. can be foliated by 
a family of space-like surfaces defined by t = x°. Let gij,i,j = 1,...,D denotes the 

3 More general form of RFDiff HL gravity was considered in Jl4|] where the action contains time and space 
derivatives of the lapse function TV according to general principles of effective field theory construction. 
However the presence of such terms has no impact on the Hamiltonian structure of given theory simply 
from the fact that the momentum conjugate to TV is not primary constraint of the theory and hence the 
Hamiltonian constraint is absent. In order to make our analysis transparent we consider the simplest version 
of RFDiff HL gravity keeping in mind that it can be easily extended to its more general versions. 
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metric on with inverse so that g%jg^ k = Sf. We further introduce the operator V,; 
that is covariant derivative defined with the metric gij. We introduce the future-pointing 
unit normal vector nf 1 to the surface £j . In ADM variables we have n° = \J — g 00 , n l = 
— g° l / \/—g 00 . We also define the lapse function N = l/y/—g 00 and the shift function 
N l = —g /g . In terms of these variables we write the components of the metric g^ u as 

500 = -N 2 + Nig tJ Nj , g 0i = N { , g^ = g tj , 



~oo = _J_ ~0i = ^ ~g _ ij N * NJ 
9 N 2 ' 9 N 2 ' 9 9 N 2 



(2.1) 



Let us now consider the general form of Hof ava-Lifshitz action 

S = ± j dtd D ^N [K^Ku - V(g)] , (2.2) 



where K%j denotes the extrinsic derivative 



K ij = ^ft&Mj ~ ViJV^ - V 3 Ni) . (2.3) 
Further the generalized De Witt metric Q l i kl is defined as 

gijkl = Itfkgjl + g U g3k) _ Xgijg kl ^ {2A) 

where A is a real constant that in case of General Relativity is equal to one. Finally V{g) is 
general function of g^ and its covariant derivative. Note also that we consider projectable 
version of HL gravity where N = N(t). 

The action fl2.2|) is invariant under foliation preserving diffeomorphism 

t'-t = f(t), x H -x l = e(t,x) ■ (2.5) 

Following (25|] we introduce U(l) transformation with parameter a(x, t) under which g%j,Ni 
and N transform as 

5 a N = , ,S agij (x,t) = , 5 a Ni(x,t) = N(t)Via(x,t) . (2.6) 



As was shown in |25| the action (2.2) is not invariant under the transformation ( |2.6] ) at 
least for Then the general procedure how to find an invariant action was formulated 



in [25]. It is based on an introducing of the scalar field v that transforms under Q2.6| ) as 

5 a v(t,x) = a(t,x) . (2.7) 
Then it turns out that the action invariant under ( |2.6[ ) can be written in the form 

S = JL J dtd D ^N((K l3 + V t V 3 u)G^ kl (K kl + V k V h v) - V(g)) (2.8) 

or in even more suggestive form by introducing 

Ni = Ni-NViU, K lJ = ^(d t g i] -V l N j -V j N i ) (2.9) 
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so that 



S = ^ J dtd D ^N [kij&Mku - V(g)] . (2.10) 
However from this analysis it is clear that v has a character of the Stiickelberg field and 



hence the symmetry (|2.6j ) is trivial. The novelty of the analysis [gjj in the formulation [26] 
is in the introduction of the additional term into action 

S v , k = J dtd D ^g( 9lj )(A - a) , (2.11) 

where 

N ■ 

a = v - N l V,v + — . (2.12) 



In the original work |2l| the function 0(g) was equal to R — Q where R is D— dimensional 
curvature and is constant. Note that in principle it is possible to consider more general 
form of Q as was suggested in p6|]. Further, a transforms under a variation as 

a'(t, x) = a(t, x) + a(t, x) - N\t, x)Vja(t, x) . 

(2.13) 

Now when we presume that A transforms under a variation as 

A'(t, x) = A(t, x) + a(t, x) - N*(t, x)Vja(i, x) (2.14) 

we immediately find that ( |2.11[ ) is invariant under a— variation. Say differently, A can be 
interpreted as the gauge field that has to be introduced when we gauge the a transformation 
P^j . More precisely, it is clear that the action ( |2.10| ) is invariant under general a(t,x) 
however as we argued this is trivial Stiickelberg extension with no impact on physical 
content of given theory. On the other hand let us presume that we want to construct 
more interesting modification of given theory when we add ( |2.11| ) without A to the original 
HL action. Now this term is invariant under a— variation on condition that a obeys the 
equation 

a(t, x) - N*(t, x)V ia (t, x) = . (2.15) 

that means that a is covariantly constant pl| and hence should be interpreted as a pa- 
rameter of a global symmetry. Gauging this symmetry means that we relax this condition 
and also introduce the gauge field A that transforms as ( 2.14| ). 



It is clear from the analysis given above that the non-relativistic covariant HL gravity 
is invariant under (2.6) for arbitrary A as was firstly stressed in [26|. Then it was argued 
that there is no scalar graviton in the perturbative spectrum about the flat background 
that makes this action very attractive since it solves the main issue of HL gravity. 

3. Hamiltonian Formalism For Non-relativistic Covariant HL Gravity 

For reader's convenience we again write non-relativistic covariant HL action 

s = ^[ dtd D ^N{k ij g i ^ kl k kl - v{ 9 )) + 

+ ^ f dtd D x^g(R)(A-a) , (3.1) 
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where we now restrict to the case when Q depends gij through the D— dimensional curva- 
ture R(gij) 4 - From ( |3.lD we find the conjugate momenta 



K 

pa ~ i Vv = — zVgG 

AC 



(3.2) 



that imply the 3 + D primary constraints 



p N ^0, p*(x)«0, $!(x) :p A (x) «0 , $ 2 (x) : p„(x) + \Vg~G(x) « . (3.3) 
Then following standard procedure we determine the Hamiltonian in the form 



H = J d D x(NH T + N l Ui + v A <5> A + v NPN + Vip 1 ] 



(3.4) 



where vn, V{, v a , A = 1,2 are Lagrange multipliers related to corresponding primary con- 
straints and where 



-Hi = -2g u V k ir kl . 

(3.5) 

Note that N and p^ do not depend on x. Now the requirement of the preservation of the 
primary constraints pjy ~ 0,pj(x) ~ 0, < I ) i(x) ~ implies following secondary ones 

d&i = {$i,H} = G = « , 

= {ptv,^} = - y rf D x^ T + \ y ^x^fvv,. « 

d D x^ T rj 
fltPi = {ft, # } = - « « . 

(3.6) 



We should stress one important issue that is related to the form of the action (3.1). At present it is 
not completely clear how the matter fields should be included in the action (3.1). Clearly we can trivially 
couple any matter field with v when we replace Ni — > Ni — NViV in all expressions containing Ni and given 
matter field. On the other hand it is an open problem how the presence of the matter field is related to 
the scalar curvature. In fact, (3.1) implies that the scalar curvature is determined solely by the function 
Q (R) = that does not depend on the matter fields. One possibility is to consider the case when Q depends 
on (R — g tJ di4>dj(f>) instead of R when the matter field is represented by the scalar field cf>. The detailed 
analysis of the action with given modification could be very interesting. 
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Now using following formulas 



{R{x),* l '(y)} = -R^(x)8(x - y) + V ! V^(x - y) - <^V fc V fe <5(x - y) 

1 — A 
XD - 1 



^g ijkl n kl - g ij v m v m g ijkl ir kl = v k (Vm kl ) + ^^r v * vv 



(3.7) 



we find that the time derivative of 3>2 is equal to 

(3.8) 



where 



(3.9) 



where generally Mij(g(x)) is a differential operator acting on tt ij that it reduces to ordinary 
multiplicative operator in case A = 1. Note that in the calculation of fl3.8|) we used following 
result 



{Pu,H} = -NV'Ui + -jViiy/gN+g) + —VV^y/gQ) « , (3.10) 



where in the final step we used the fact that the result is proportional to the constraints 
Hi and ~ 0. In the same way we find that 



2 {^'/ rf '° xiV ' I/V < v i 7r< ' 7 '( x )) 
J d D xNVivH\ yfg~g\ ~ y/g^gsTv ps . 



(3.11) 

Let us review constraints that we derived at this stage. We have following set of 
secondary constraints ~ , $2 ~ ~ and one global T = J d xHt ~ 0. Note 
also that = $1 — <&{ ~ that according to ( |3.10| ) is the first class constraint. Then the 
total Hamiltonian takes the form 

H T = J d D x{NU T + N l U t + v A PA + v NPN + v iP l + v u Pu + vjM 1 + Wjj*") . 

(3.12) 

where vn ,Vi,v"^,Vjj,Vjj are corresponding Lagrange multipliers. Note that we included 
the expression [A — N l Viv + y V*Vj^) into definition of the Lagrange multiplier vjj. 

As the final step we analyze the stability of the secondary constraints. Let us begin with 
the constraint %i . It is convenient to extend these constraints by appropriate combinations 
of additional constraints p u ps 0,p A ps so that 

Hi = -2g lk ViTr kl + diApA + d t up v . (3.13) 



-7- 



Then Ts(iV*) = J d D is.N i 'Hi is generator of the spatial diffeomorphism that is clearly 
preserved during the time evolution of the system since the Hamiltonian is invariant under 
spatial diffeomorphism. Further, p u ~ is preserved during the time evolution of the system 
according to fl3,10| ). On the other hand the time evolution of the constraint ~ is 
equal to 

dM 1 = {vl'tHr} « J d D x (^^ / (x) + «? / (x){*P,^(x)} 
^f / (x){${ / ,^(x)} = 0. 



(3.14) 



Since 



fcS'br)} = A4y(y) ({<&{' (x),^(y)}) 



(3.15) 



we find that the equation ( 3.14| ) implies that vjj = 0. In the same way the requirement of 
the preservation of the constraint implies 

d t $2 = N{*2,T} + J d D x^ / (x){$ 2 ,$ 1 (x)} = 

(3.16) 

that due to the fact that {<E>2,T} ^ and ( |3.15| ) allows to determine vjj as a function of 
the canonical variables. In other words, ${ J and $^ are the second class constraints 5 

We see that the requirement of the preservation of the secondary constraints does not 
imply additional constraints so that we obtained following constraint structure. We have 
first class constraints Hi « 0,p v « 0,pi « 0,p^ « together with two global first class 
constraints pn ~ 0, T ~ 0. Then we have two second class constraints . The 

detailed discussion of these constraints will be given in the next section. 

In this section we performed the Hamiltonian analysis of non-relativistic covariant 



HL gravity in the formulation presented in [26] and we showed that it leads to the same 



5 The nontrivial property of given theory is that the Poisson bracket between the second class constraints 
depend on the phase space variables so that it is possible that it vanishes on some subspace of phase space. 
In order to analyze this issue we should explicitly determine the form of this Poisson bracket and after 
some algebra we find {$"(x), $^(y)} = A(R, Rij, x, y) + (1 ~^~ D) V a V*Vj V J £(x - y)- We see that due 
to the second therm this Poisson bracket is non-zero on the whole phase space on condition when A / 1. 
On the other hand in case when A = 1 we find that this Poisson bracket vanishes for the subspace of the 
phase space where Rij =0. However when Rij — we see that $2 is preserved during the time evolution 
of the system and hence it is not necessary to impose additional constraint $2 7 ~ 0. Then &1 1 w is the 
first class constraint. The gauge fixing of given constraints implies that all metric components and their 
conjugate momenta are non-propagating degrees of freedom and hence the theory on the subspace Rij = 
is topological with no local degrees of freedom. 
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structure of the constraints as in the original proposal [25|. Note that our analysis is valid 



for general A with agreement with [26|. Further, as was shown in [25] the number of physical 
degrees of freedom is the same as in the General Relativity even if the constraint structures 
of these two theories are different. On the other hand the non-relativistic covariant HL 
gravity has an additional global Hamiltonian constraint. However when we consider RFDiff 
invariant HL gravity as the starting point for U(l) extension of HL Gravity we find theory 
with the same content of physical degrees of freedom as in non-relativistic covariant HL 
gravity with additional important difference which is an absence of the global Hamiltonian 
constraint. 

4. Non-Relativistic Covariant RFDiff HL Gravity 

RFDiff invariant Hofava-Lifshitz gravity was introduced in [14] and further studied in 



1 23]. This is the version of the Hofava-Lifshitz gravity that is not invariant under foliation 



preserving diffeomorphism but only under reduced set of diffeomorphism 

t' = t + St , 5t = const , x H = x i + f (£, x) (4.1) 

As was argued in [23] the simplest form of RFDiff invariant Hofava-Lifshitz gravity takes 
the form 

S = ^ J dtd D ^{k i3 G^ kl k kl - V(g)) , (4.2) 

where 

Kij = \(d tgij - ViNj - VjNi) . (4.3) 



Note that this action differs from HL gravity action (2^) by absence of the lapse N and 
by replacement of the extrinsic curvature with given above. This action is in- 
variant under RFDiff symmetries ( |4.lD that is reduced with respect to foliation preserving 
diffeomorphism . 

In order to find the U(l) extension of given theory we introduce the field v and replace 
Ni with Ni as 

N i = Ni- Viu . (4.4) 
Then it is again easy to see that the action is invariant under transformation 

Nl(t, x) = Ni(t, x) + Via(t, x) , u'(t, x) = u(t, x) + a(t, x) . (4.5) 

Clearly this replacement is as trivial as the one performed in the projectable version of 
Hofava-Lifshitz gravity. Then following the same procedure as in section (|2[) we find the 
action in the form 

Srfd = ^ f dtd D ^(K l3 g^ kl K kl - V{g) + G{R){A - a)) , (4.6) 



where 



Kij = -(d t gij - ViN 3 - - VjNi + ViVjV + VjViu) . (4.7) 



-9- 



Clearly this action is invariant under (|4.1| ) and under (4.5). Further A and a transform as 
scalar under ( [4.1| ) 

A'(t', x') = A(t, x) , a'(t', x') = a(t, x) (4.8) 



Note that the action Q4,6| ) can be derived from non-relativistic covariant HL action by 
setting N = 1 and hence one can expect that these theories describe the same local physics. 
However the my difference between these two formulations emerges when we perform the 
Hamiltonian analysis of the action (|4.6| ) . 

As in previous section we find the primary constraints 

Pi (x) « , $! : p^(x) « , $ 2 : p„(x) + -^^(x) (4.9) 

and the relation between K{j and conjugate momenta 7t lJ 

Kij = ~jf^ kl ■ ( 4 - 10 ) 
Then it is easy to find the total Hamiltonian in the form 

H = j d D *{U T + N l Hi + v% + v NPN + vd) - 



i J d D *Jgg{R){A- N^iu +\v i V i v) , 



where 



\J9 K K 



-Hi = -2g u V k 7T kl 



d m = { Pi , H} = -Ui - $f « -Tfc « . 



(4.11) 



(4.12) 



The requirement of the preservation of the primary constraints Pi(x) ~ 0, 3>i(x) ~ implies 
following secondary ones 



(4.13) 



In case of the preservation of the constraint $2 we proceed as in previous section and we 
find 

(4.14) 
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where 



&1 1 = -R 



A „ 1 - A 



(4.15) 



and where generally Mij(g(x)) is a differential operator acting on tt^ that it reduces to 
ordinary multiplicative operator in case A = 1. Note also that p v ~ is the first class 
constraint. 

Following general analysis of constraints systems we introduce the total Hamiltonian 
in the form 



H T = J d D x(n T + rfHi + v A p A + v u Pu + Vip 1 + v}M 1 + vjM 1 ) 



(4.16) 



As the final step we should perform the analysis of the secondary constraints. However 
this was done in previous section so that we do not repeat it here. 

Let us now discuss the second class constraints 3>{ J , $| . According to standard analy- 
sis these constraints have to vanish strongly and allow to solve for two phase space variables 
as a functions of remaining physical phase space variables that span the reduced phase 
space. However solving these constraints in full generality is very difficult. On the other 
hand it is easy to see that in linearized approximation these constraints can be solved as 
h = , 7r = where h is the trace part of the metric fluctuation and tt is its conjugate 
momenta. 

Even if we cannot solve these constraints explicitly in general case we can still determine 
the number of physical degrees of freedom. To do this note that there are D(D + 1) gravity 
phase space variables gij,^, 2D variables Ni,p % , 2 variables A,p A and 2 variables v,p v . 
In summary the total number of degrees of freedom is No.o.f = D 2 + 3-D + 4. On the 
other hand we have D first class constraints Hi ~ 0, D first class constraints pi ~ 0, 2 first 
class constraints p v ~,pa ~ and two second class constraints ^l 1 ,^^ 1 - Then we have 
Nf.c.c = 2-D + 2 first class constraints and N s . c , c . = 2 second class constraints. Then the 
number of physical degrees of freedom is [plf| 

No.o.f. ~ 2JV/.CC - iW = D 2 — D — 2 (4.17) 

that exactly corresponds to the number of the phase space physical degrees of freedom of 



D + 1 dimensional gravity. For example for D = 3 the equation (4.17) is equal to 4 which 
is the number of phase space degrees of freedom of massless graviton. 

In summary the Hamiltonian of non-relativistic general covariant RFDiff HL gravity 
gives the appropriate number of physical degrees of freedom of gravitational theory without 
introducing global Hamiltonian constraint. There is also another interesting aspect of given 
theory which is its non-trivial symplectic structure. In fact, let us denote the constraints 
3>i2 as wriere A,B = I, II so that the Poisson bracket between constraints can be 
written as 

{^(x),$i/(y)}=A AB (x,y). (4.18) 
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From the structure of these constraints we find that the matrix Aab has following structure 



A AB (x,y) = f ° * | , (4.19) 



* * 



where * means non-zero elements. Then the inverse matrix (A 1 )ab has the form 

(A-y 13 ( m ) . (4.20) 

Now we observe that 




0fc(*),*i J (y)} = ,{ 5ii (x),$f(y)} ^0 , 
{^'(x), $^(y)} + , {^(x), (y)} . 



(4.21) 



Then we find that the Dirac brackets between canonical variables take the form 



{9ij(x),9kl(y)} D = -J dzdz / {^(x),^(z)}(A- 1 )^(z,z'){^ / (z'),^(y)} = 0, 

-J dzdz'{^(x),^(z)}(A- 1 )^(z,z'){^ / (z / ),vr H (y)} = fF H (x,y) , 
{^(x),vr fc '(y)} D = {^(x),vr fc '(y)}- 

| dzdz'{^(x),^(z)}(A- 1 )^(z,zO{^ / (z'),^(y)}=^(x,y) , 

(4.22) 



where the matrix O depends on phase-space variables according to ( f4.20j ) and (4.21). Hence 



the non-relativistic covariant RFDiff HL gravity has well defined Hamiltonian formulation 
with symplectic structure that generally depends on phase space variables. Note however 
that in case of the linearized approximation one can choose the constraints in such a 
way that the Dirac bracket coincides with the Poisson bracket. Explicitly, in linearized 
approximation the second class constraints can be chosen as h = 0, ir = as follows from 
the analysis given above. These constraints have vanishing Poisson brackets with remaining 
dynamical variables and consequently the Dirac brackets between physical phase space 
variables coincide with Poisson brackets. 

As the final remark we again emphasize the important point that the Hamiltonian of 
non-relativistic covariant RFDiff HL gravity does not vanish on constraint surface. This 
is the similar situation as in case of the Hamiltonian of the healthy extended HL gravity 



1 17, [lg] which is however in sharp contrast with the Hamiltonian of General Relativity. As 
we argued in these papers this fact has a strong impact on the definition of observables in 
healthy extended Hof ava- Lifshitz gravity or in any theory of gravity where the Hamiltonian 
is not given as linear combination of constraints. Since the discussion presented in [|l7], 18] 
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can be applied in case of the non-relativistic covariant RFDiff HL gravity as well we are 
not going to repeat it here. Instead we recommend these papers to reader that is interested 
in these problems. 

It is important to stress that the fact that the Hamiltonian is not given as linear 
combination of constraints has an important consequence for the stability of given theory. 
Explicitly, it is well known that some massive gravities are unstable since the Hamiltonian is 
not bounded from bellow. Alternatively, the instability of given theory is also indicated by 
presence of the ghosts (fields with wrong sign of kinetic term) in the fluctuation spectrum. 
However we believe that this is not the case of non-relativistic covariant RFDiff HL gravity 
even if the full analysis has not been done yet. The crucial fact is the absence of the scalar 
graviton in the fluctuation spectrum that implies that RFDiff invariant HL gravity has 
well defined Hamiltonian that is positive definite at least in linearized approximation. In 
order to understand the properties of the Hamiltonian of non-linear theory we should solve 
the second class constraints and express the Hamiltonian in terms of physical modes only. 
However as we argued above this is very difficult task and hence the problem of the stability 
of the general non-relativistic covariant RFDiff HL gravity has not been fully addressed. 
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